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Abstract 

We extend the recent gravitational decoherence analysis of Pikovski et al. to 
an individnal mesoscopic system with internal state characterized by a coher¬ 
ent snperposition of energy eigenstates. We express the Pikovski et al. effect 
directly in terms of the energy variance, and show that the interferometric 
visibility is bonnded from below. Hence nnlike collisional decoherence, the 
visibility does not approach zero at large times, althongh for a large system 
it can become very small. 


In a recent interesting article, Pikovski et al. [ij show that there is a nni- 
versal decoherence-like effect for systems in a varying gravitational potential. 
They focns on a system consisting of a large nnmber N of thermally excited 
harmonic oscillators, and calcnlate a formnla, with which we agree, for the 
time dependence of the interferometric visibility. Their formnla depends on 
the snbsystem nnmber and temperatnre T as N^^'^ksT (with ks the Boltz¬ 
mann constant), and as they note this indicates that the internal energy 
variance is the relevant system attribnte for their effect, which is explicit 
in their earlier paper discnssing gravitational decoherence in a two state 
system. Two recent papers j^, have noted that the Pikovski et al. effect 
vanishes in a freely falling frame, and so it fnrnishes yet another interesting 
example of the snbtle interplay of qnantnm theory and gravitational physics. 
In this paper, however, we conhne onrselves to the frame of an asymptoti¬ 
cally (with respect to the Earth) inertial observer; here the nonrelativistic 
Newtonian limit snffices to give correct answers. 
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In this note we extend the derivation of [l(] to a single mesoscopic system 
characterized by a general coherent snperposition of energy eigenstates, and 
show directly that the gravitational reduction in interferometric visibility 
depends, for small times t, on the energy variance AE. We also show that 
the visibility is strictly bounded from below, and does not approach zero for 
large times. When extended to a large collection of N independent meso¬ 
scopic subsystems, with thermal energy variance AE our result 

maps to that of Pikovski et ah, and the lower bound vanishes exponentially 
for large N. However, this exponential vanishing is critically dependent on 
the independence assumption, and may not apply to strongly coupled sub¬ 
systems. 

To keep this article self-contained, we present a complete derivation, start¬ 
ing from the underlying gravitational physics. According to general relativ¬ 
ity, gravity couples to matter via the stress-energy tensor T^j,(x); in the non- 
relativistic limit and to first order in weak gravitational fields, the interaction 
Hamiltonian is 

Hc = ^ J d^a;0(x)Too(x) , (1) 

where </>(x) is the Newtonian potential. If the energy density is well localized 
in space with respect to the distances over which 0(x) varies appreciably, we 
can bring (/>(x) outside the integral to obtain Hq = where x is 

the mean position of the matter content and H is the total Hamiltonian of 
the matter distribution. In the case of the Earth’s gravitational field, where 
4>{x) = gx, with g the Earth’s gravitational acceleration and x the vertical 
distance above Earth’s surface, we therefore have 

Ho = gxHlc^ . (2) 


The total H is given by 


H 


me + H 


center-of-mass kinetic 




( 3 ) 


with Hint the system internal energy. The dominant term is the rest mass mc^, 
which produces the standard Newtonian term V = mgx in the Schrodinger 
equation, the quantum effects of which have been measured in neutron in¬ 
terferometry 1^. The internal energy contribution in Eq. (|3]) introduces 
new effects, specifically a phase shift in the time evolution of the internal 
energy states, which is the origin for the time dilation decoherence discussed 
by Pikovski et ah In the calculation that follows, we neglect the mc^ term 


2 


in the energy, as well as the kinetic energy associated with center of mass 
motion, since these do not contribute to the internal energy variance and so 
drop out of the final formula. 

Let us now consider a single system and calculate for this the analog of the 
Pikovski et ah effect. We assume that the total initial state of the system 
is |'0(a;,O)) = \'i1^cm{.x)) 0 iV'mt), where \'i1^cm{x)) is the center-of-mass wave 
function, which is assumed to be well-localized in space around the position x, 
while It/'int) is the internal dynamics, which we decompose as a superposition 
of eigenstates \n) of the internal Hamiltonian: The state 

at time t is thus 


\ih(T. 6d I'll), (t. 



(4) 


We see that due to the presence of the Earth’s gravitational held, the internal 
states acquire a phase which depends on the position of the center of mass 
with respect to the Earth. For a single system in an energy eigenstate, 
the gravitational effect appears only as an overall phase of the internal wave 
function and does not change its magnitude; the more interesting case is when 
two or more energy eigenstates appear in the sum over n in Eq. (jl]), so that 
the gravitational effect cannot be factored out as an overall phase. Note that 
even in this more complex case, the gravitational effect is reversible, in that if 
the system is sent through a gravitational field, and then through an equal but 
reversed gravitational field, the effect is eliminated. This already indicates a 
contrast with standard collisional decoherence induced by interaction with a 
chaotic environment, which cannot for all practical purposes be reversed. 

We now consider a system that is in a superposition of two different posi¬ 
tions in space, |'0 cm( 2^)) —t [|'^cm(2:i)) + |'^cm( 2 : 2 ))]/\/ 2 , with the two center- 
of-mass states well localized with respect to Ax = Xi — X 2 , and practically 
orthogonal. Then the internal states will entangle with the center of mass 
and the reduced density matrix ff^ of the center of mass, obtained by tracing 
the density matrix over the internal degrees of freedom, will dephase. More 
precisely, the off-diagonal element = (' 0 cM(a^i)|p‘^'“|V’cM(a^ 2 )) evolves in 
time as follows. 




n 


We will be interested in calculating the interferometric visibility, which here 
is twice the absolute value |Pi|^(t)|. This is unchanged when we multiply 
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Pifit) by any phase factor, and it is convenient to choose this phase factor as 
follows. Let us dehne E = \cn\^En as the mean internal energy, and dehne 
= En-E, so that Y.n IcnpAEn = 0 and (AE)'^ = Y,n |cnP(AE„)^ gives 
the square of the internal energy variance AE. Then we get the same absolute 
value \Pi^{t)\ if we replace pi^{t) by pi^{t) = p1f{t). 

Expanding in powers of f, we get 


2pf2“(*) = E 


AEnAx 
ig——t 


hE 


1 ^{AEn)^Ax^ 
2^ 


+ 0{t^) 


= i-e/ti + o{e), 


( 6 ) 


where we have introduced the phase evolution time scale (which governs the 
Pikovski et ah effect). 


y/2hc^ 

gAE\Ax 


(7) 


The phase evolution time depends on the spatial separation of the two cen¬ 
ter of mass components as well as on the internal energy variance, and the 
behavior of Eq. ([S]) depends strongly on the system complexity. For a single 
system consisting of a superposition of only two energy eigenstates, 2|p°^(t)| 
oscillates in time, on a time scale of order of magnitude t£). For a single sys¬ 
tem consisting of a superposition of many eigenstates, over large time scales 
the behavior of the visibility is still oscillatory, and is governed by the lower 
bound derived below. But for a complex system the behavior of the interfer¬ 
ometric visibility for small times f is of interest in understanding the 

Pikovski et ah effect, and is given by 


2\p-f{t)\c^i-e/ti 


( 8 ) 


and for a system composed of N identical subsystems, the visibility for small 
times is given by 

‘^\Pi2ii) \ ^ (1 - ^ exp(-iVtV^D) • (9) 


This can be rewritten as 


2|p-2"Wl^exp(-tV4j 

with 

'/2hc^ 

” VNgAE\Ax\' 


( 10 ) 

( 11 ) 
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in which the effective energy variance Ai? is multiplied by \/N. Note that the 
exponentially vanishing behavior of Eq. ([2]) for large N is critically dependent 
on the assumption of N independent systems. For a strongly coupled N 
particle system with A.Ej\f ~ N^^'^AEi, without further input we could only 
conclude that the right hand side of Eq. ([9]) is replaced by 1 — Nt^/t\j^ 
for small t (i.e. Eq. ([2]) with AE replaced by AEj^, the total internal 
energy variance), which does not imply exponentially vanishing behavior of 
the visibility for large N. A similar cautionary remark applies to the lower 
bound on the visibility of a system composed of N independent subsystems 
derived in Eq. ffT7|) below. 

If the initial internal state is not a pure state, but a statistical mixture of 
states with probabilities Pa, then Eq. ([7]) still holds, with; 

( A £)2 = ( 12 ) 

a^n 

which measures internal energy fluctuations both of ‘quantum’ origin due 
to the eigenstate superposition with amplitudes c„, as well as of ‘classical’ 
origin due to the statistical probabilities Pa- Pikovski et al. consider the 
case of only thermal fluctuations, in which case AE = y/NksT, where N 
is the number of degrees of freedom, ks Boltzmann’s constant and T the 
temperature, and one recovers Eq. (3) of their paper. This also agrees with 
Eq. fITT]) above for the case of N independent subsystems with individual 
energy variance AE = ksT. 

We now show that the effect produced by the coupling of the internal 
degrees of freedom to the gravitational field, being a sum of phase shifts, 
does not correspond to decoherence in the usual sense. Forming the absolute 
value squared of the visibility, we have 

{2\p^f{t)\Y = (13) 

n,m 

which implies that the average of the visibility over a large time interval is 

rii+oo ^ > O’ 

°° Jo ^ 

where we have assumed AEn 7 ^ AE^ for n ^ m. When there are energy 
degeneracies, this formula still holds, with |cnp the sum of the absolute value 
squared coefficients of all states with the same energy. 
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Since 'Yhn ~ lower bound of Eq. flT^ is always a number 

smaller than 1, but greater than 0. A stronger statement can be made when 
the number of states in the superposition is a finite number L. Expanding 
out the inequality 

L L 

m=l n=l 

and using state vector normalization J2n=i ~ ^ counting 

ELi 1 = ^ gives 

L 

'Y^\cn\^>L~^ ■ (16) 

n=l 

Equations flT^ and flT^ mean that |pi“(t)| cannot vanish for large times. 
This is in contrast to usual collisional decoherence, where the effect of a 
single interaction is of the form —)> Ci 2 Pi^ with |ci 2 | < 1, and so after Ni 

interactions, the off-diagonal element approaches zero as |ci 2 |^b When W is 
linear in time this gives an exponential vanishing of the visibility as a function 
of time. By way of contrast, the gravitational decoherence calculated above 
gives a visibility that vanishes exponentially as the number of independent 
subsystems N approaches infinity, but which has a Gaussian decrease in time 
only for small times. For N independent subsystems, the lower bound of Eq. 
(ITT)) is replaced by 

El'"!")" . (V) 

n 

which since 'Yhn \^ri\^ ^ 1 approaches zero exponentially as the subsystem 
number N approaches infinity. 

As a concrete illustration of the smallness of the lower bound for large 
systems, consider a cube 10“^ cm on a side, containing roughly 1000 atoms, 
with of order 1000 acoustical vibration modes. Since these modes are ap¬ 
proximately independent, we can apply Eq. 03 ; assuming that the cube 
state is prepared so that the average mode is in a superposition of 3 states, 
we can also use Eq. (IT^ with L = 3 as an estimate for Yhn giving for 
the lower bound on the long time average of the visibility 

lim 1 f (2\p<S(t)\fdt > 3-™", (18) 

T^+oo 1 Jo 

which is zero for all practical purposes. With this same model in mind, it is 
also instructive to consider the competition between gravitational decoher- 
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ence and standard collisional decoherence. For collisional decoherence, in the 
limit of small spatial snperpositions, the analog of Eq. flTU]) is Q 

= |pf2^(0)|exp(-t/fcoii) , (19) 


with 


^Coll 

A 


1 

A|Aa:P 

na{q^v)AV 

3A2 


( 20 ) 


In this formula, n is the density of scatterers that scatter from the decoher¬ 
ing object with cross section a, and q and v are respectively the scatterer 
momentum and velocity. If the scattering particles have mass m and form a 
thermal bath at temperature T, then 

(g^'y)Av = 4(m/7r)^(2fcBT)t . (21) 


Assuming that the decohering object is in thermal equilibrium with the bath, 
we can take its temperature also as T, and so in the formula for tN^ ^6 can 
take AE = y/NksT. Combining the various formulas, the condition for 
tNo < ^Coii can be written as 


3(A7r)2 hg 

16 c^\Ax\a{mkBT)^ 


( 22 ) 


Taking as an example the 10 atom cube and assuming the scatterers are 
nitrogen molcules at room temperature, we have as inputs for a numerical 
estimate of decoherence of a superposition in which the cube center is dis¬ 
placed by the cube diameter, 

a =10-^W2 , 

N =1000 , 

*.BT=ieV . 

|Aa;| =10 '^cm , 
m =14 X 10®eV/c^ , 

g =981cm/s^ , 

A =6.6 X 10“^®eVs , 

c =3 X 10^°cm/s , 

(23) 
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the inequality of Eq. fl2^ becomes 

n < 1.2 X 10“^cm-3 . (24) 

This is a density of 10“^“^ of atmospheric density at standard temperature 
and pressure, and corresponds to a vacuum presently unattainable in the 
laboratory. So under normal laboratory conditions, collisional decoherence 
occurs on a more rapid time scale than gravitational decoherence. 

To summarize, we have given a ^neralized and simplified derivation of the 
gravitational decoherence effect of [ij] that applies to a single mesoscopic sys¬ 
tem containing an arbitrary superposition of energy states, and have shown 
that there are important limitations on the size and applicability of the effect. 
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